Stimulated Raman scattering of water maser lines in astrophysi- 
cal plasmas 



R. T. Gangadhara 

Indian Institute of Astrophysics, BangaIore-560034, India 
Shuji Deguchi 

Nobeyama Radio Observatory, National Astronomical Observatory, Minamimaki, Minamisaku, 
Nagano 384-13, Japan 

H. Lesch 

Institut fiir Astronomie und Astrophysik der Universitat Miinchen, Scheinerstra/3e 1, 81679 
Miinchen, Germany 

Radiative transfer equations are derived and solved for the stimulated Raman scattering of 
water maser lines in the astrophysical plasmas with electron density of about 10^-10^ cm^'^. 
In stimulated Raman scattering, the energy of water maser line is transferred to the side band 
modes: Stokes mode and anti-Stokes mode. The Stokes mode is easily produced by backward 
Raman scattering while the anti-Stokes mode is created by the interacting intersecting masers 
in the plasma. The intensity of the Stokes mode is higher than that of the anti-Stokes mode. 
These side band modes are proposed as explanation for the extreme velocity features observed 
in the galaxy NGC 4258. The threshold value of the brightness temperature for the Raman 
scattering is about 10"*^^ — 10^^ K, and it is satisfied in the case of NGC 4258. 
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I. INTRODUCTION 

Stimulated scattering processes become important in astrophysical plasmas when the incident 
radiation on plasma has a very high brightness temperature (Tf, > 10^^ K) ^'2'^'^. In practice, 
line radiation of water masers in the galactic nuclei and continuum radiation from quasars are 
considered to have such an high brightness. Stimulated Raman and Compton scatterings will 
occur in the high and low density plasmas, respectively, depending upon whether the collective 
effect of electrons works or not ^. Though a number of studies on the stimulated Raman scattering 
have been made, the transfer of line radiation with shifted frequency is not well understood. 
Fernandez and Reinnisch^ have made an attempt made to formulate the theory of stimulated 
Raman scattering with radiative transfer equation for specific intensities. 

In this paper, we have made an attempt study stimulated Raman scattering of water maser 
line with shifted frequency. A practical application of this theory is aimed to explain the extreme 
high-velocity features of water maser lines at 22.235 GHz in the galactic nucleus of NGC 4258. 
The high-velocity features of this line appears separated by 1000 km s""^ from the main com- 
ponent of the water maser line. Though the recent Very Long Baseline Interferometry (VLBI) 
observations showed some evidence of this component arising from the rotating disk, all the 
observations (for example, the weakness of the components) are not necessarily well explained 
by the rotating disk model ^. The beaming of Raman maser is discussed in Sec. II. The radiative 
transfer equations for stimulated Raman scattering are derived and solved in Sec. III. Further 
we discuss the limitation and applicability of our model to the astrophysical maser phenomenon 
in Sec. IV. 

II. BEAMING OF RAMAN MASERS 

The principal mechanism behind the beaming of astrophysical masers is the frequency change of 
the resonance radiation. As usual for astrophysical masers, the Raman masers are also expected 
to be strongly beamed. The masers are amplified in the direction with smallest gradient in the 
resonance frequency where the optical depth of the resonance radiation becomes maximum. For 
the case of Raman masers, the change in the resonance frequency is produced by two effects: 
change of plasma frequency due to temperature and density variations, and Doppler effect due 
to the gas flow in plasma. 

The frequencies of incident and scattered waves are related by 

UJi=Uo-i03, (1) 
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where cOg, loi and co^ are the frequencies of incident, scattered and longitudinal waves, respec- 
tively. The frequency of the longitudinal wave is determined by the dispersion relation 

iol = iol^ + 3v^kl (2) 

where = ksT /me and lJ^^ = AirNeC^ /nie, e and rUe are the charge and mass of electron, Ne 
and T are the density and temperature of electron plasma, and is the Boltzmann constant. 

For the case 0^3 <C oJo, we can approximate ks = 2u>ocos0s/c and ^3 = (tt — 9i)/2, where 83 
is the angle between k^ and ^3, and 9i is the angle between kg and ki, and c is the speed of 
light. Therefore, the dispersion relation can be transformed as 

In astrophysical plasmas, density and temperature of the plasma are not constant in space, 
and hence the resonance frequency varies along the propagation path of the waves. 

In addition, the Doppler effect due to the gas flow in plasma changes the frequency, which 

can be expressed in the rest frame as 

uji = [iOo{l + v.ho/c) - UJ3]{1 - v.fii/c), (4) 

where v is the velocity of the plasma in the rest frame, and ho and hi are the unit vectors 
representing the propagation directions of incident and scattered waves, respectively. 

The frequency shift along the co-ordinate s, which is taken in the direction of scattered wave, 
can be expressed as 

^ = Uo{hi.V){v.ho/c) - wio(ni.V)('U.ni/c) - (ni.V)w3, (5) 

where ujiq is the original value of frequency u}\. The first two terms on the right-hand-side 
represent the frequency change due to Doppler motion of the gas, and third term is due to the 
density and temperature changes in plasma. 

The transfer of line radiation in the moving medium has been well studied and an approxi- 
mation of large velocity gradient has been developed. For example, for a spherically expanding 
cloud, the Doppler term (in this case u}o = ^i) can be written as 

(^) (6) 

V as /Doppler c dr 

where wc have taken the z axis parallel to the direction of incident radiation, which is considered 
to be the direction of maximum optical depth {9i = 0) . The frequency of the scattered radiation 
is not influenced by the Doppler motion of the matter in the case of pure-forward scattering. 



The density-temperature gradient term can be written in unit of (lOq/c) (dv/dr) as 



= — ^[cos 9i{az + bz{l — cos^i)} + sin^i cos(f)i{ax + bx{l — cos^i)} 
ds / c dv 

+ sin^i sin^ijoj, + 6^(1 — cos^i)}], (7) 

where vectors a and b are the non-dimensional parameters, which are proportional to the density 
(a) and temperature (b) gradients. 

III. STIMULATED RAMAN SCATTERING OF MASER 

The fundamental set of equations for the study of stimulated Raman scattering in the plasma 
medium are 

^ + V.iN,v) = 0, (8) 

_ + {T.Vv = [E + -VXB] - —vj^VNe - i^cv, 9) 

ot nie \ c J Ng 

V.E = 47re(iVi - N^), (10) 
Vx^.-if, (U) 
-J 47r^ IdE 

The Eqs. (8) and (9) are the electron continuity and momentum equations, respectively. The 
Eq. (10) is the Poisson equation. The last two Eqs. (11) and (12) are Maxwell's equations. Here 
E and B arc the electric and magnetic fields of the waves involved in the three-wave process. 
Taking curl on the both sides of Eq. (11) and substituting Eq. (12) for V x 5, we get 

^ + c2Vx(Vx^)-h47r^ = 0. (13) 

Let Ne = No + dN and Ni = No be the electron and ion densities, respectively, where SN is the 
electron density perturbation induced by the ponderomotive force of the radiation field, and Np 
is the stationary background neutralizing ion density. 
The current density is defined as 

j = e{NiVi-NeVe). (14) 

Taking the time derivative of j and substituting into Eq. (13), we get 

d'^E 2^ , ^dv ^ d(6Nv) ^ 

^ + c'V X (V X £;) - 47reiVo^ - 47re^^ = 0. (15) 
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Prom Eq. (10), we find 

SN = —V.E. (16) 

47re 

Using Eqs. (9) and (16), Eq. (15) can be written as 

a2 Tp AT 

+ c^V X (V X ^) + ujI^E - 3-^4V(V.^) = -47reNo{v.Vv + iycv) 



e 



^le,^ ^^ , d{6Nv) 

^{vxB) + 4:Tre \' . (17) 

c ot 

Now, using vCV.v) = (1/2) Vv^ —vx {V xv) and V xv = eB/{mec), ^, Eq. (17) can be written 
as 

__ _ c2 - - ^4 V(V.£;) + = - l-KeNoVv" - A-kcNoVcV. (18) 

We use plane-wave approximation ^ such as 

E{r, i) = I E t)e'^'''*-^'-^ + E*(r, t)e-^('^'*-^-*^} (19) 

j=0-3 

and a similar approximation for the electron motions: 

^(f^ i) = ^ E i)e^('^'*-^-'^ + <(f, i)e-^('^^*-^-^^}. (20) 

i=0-3 

In Eqs. (19) and (20) the amplitudes Ei and Vi are the slowly varying functions of space and 

time, where i = 0, 1, 2 and 3 stand for incident (pump) wave (fcojCJo), Stokes mode {ki,u)i), 
anti-Stokes mode (^2,1^2) and longitudinal plasma wave (fcsjtJa), respectively. 

The stimulated Raman scattering instability excites resonantly when the following phase 
matching conditions are satisfied: 

wi=Wo-W3, ki = ko-k3, (21a) 

U2 = 0J0 + ^3, k2 = ko + ks- (21b) 
The growth rate of instability is given by ^^'^^ 

7o = ^f^V'', (22) 



2 \uji 

where Vg = eEg/ (mgCJo) is the quiver velocity of the electrons in the incident radiation field. 

Using the first-order approximation ^ to the electron motion we find the transfer equations 
for electric field amplitudes: for the incident radiation (pump) 

QE dE 

"df ^ ""d^ " " "^^1^3 cos - E2E*s cos V2), (23) 



for the Stokes mode 



for the anti- Stokes mode 



"df ^ "'d^ " ^ "^"^3 cos Vi, (24) 



'W ^ ""a^ " ""^'^^ " cos V2 (25) 

and for the longitudinal wave 

dE-^ dE% 

-Qf+''9-g^ = -C3^3 + MEoE^ COS ^1 + E2E; COS ^2), (26) 

where q = {up^/uJi)^Vc (i = 0, 1, 2, 3) represents the free-free damping process, and the electron- 
ion collision frequency Vc = 47rnee'^ In A/mgvf . The coefficients arc a = k^e/{AmeiOo), Q^s = 
iOp^a/luJoijJs)- Here Vg = Svj^ks/uJs is the group velocity of longitudinal wave and tpi {i = 1 or 2) 
is the angle between Eg and Ei or E2. The co-ordinates Xi (i = 0, 1, 2, 3) represent the direction 
of propagation of each wave. 

Only the frequency downshifted component (Stokes mode) is excited, and the upshifted 
component (anti-Stokes) is not amplified at the general scattering angles even if it is resonant 
It should be noted that the second terms on the right-hand-side of Eqs. (24) and (25) are 
plus and minus. This means that only the downshifted waves grow, but the upshifted waves never 
grow. However the resonance can occur for either of the upshifted or downshifted components 
in general, but both the components appear only in the case of pure forward scattering. 

For the longitudinal wave. Landau damping is more efficient than the binary coUisional 
damping. The Landau damping rate of the plasma wave is given by 

In the steady state, Eqs. (23), (24) and (26) reduces to 

9^ = _2^E^_-E,E,, (28) 
0X0 c c 

^ = -^Eo + -EoEl 29 
0x1 c c 

Es = -^EoEl, (30) 

r + C3 

where we have assumed incident wave and scattered wave are polarized in the same direction 
(V'l = 0). Eliminating E3 between the Eqs. (28) and (30), we find 

P' = --Eo-^^,\E,\^E^. (31) 

0X0 C C{r + C3) 



The complex conjugate of Eq. (31) is given by 



-K - -7^^ El. (32) 



dXo C ° c(r + C3 

Multiplying Eq. (31) by and Eq. (32) by Eg, and adding, we get 



dio 2co J 2Q:a3 

-w- = lo - /p , J ilo- (33a) 

OXo C C(r + C3) 



Similarly from Eq. (29), we find 

dh 2ci 2aa3 

^ = h + -7p- (336) 

0x1 c c{T + C3) 

where Iq =\Eo\'^ and Ii =|£^ip are the intensities of the incident and scattered waves. In 
Eq. (33) the first term on the right hand side represent the free-free collisional damping. While 
the second term in Eq. (33a) represents stimulated absorption of the incident radiation and in 
Eq. (33b) it represents the stimulated emission of the scattered radiation. 

A. Solution of equation (33) 

When the amplitudes of and Ii are not too large, we solve Eq. (33) numerically. But when 
If, and Ii are too large we can ignore the collisional damping terms and find the analytical 
solutions. 

The direction of propagation of the incident and the scattered waves in the thick and the 
thin plasma media are shown in Fig. 1. Let be the angle between the coordinates Xg and xi, 
then we have Xp = xi cos 9, and Eq. (33) becomes 



—— = — cos(y) 

ds 



(34a) 



where ds = Agdxi, Ag = 2co/c, Bg = aa3/{c{T + C3)}, Ai = 2ci/c, and Bi = Bg. 

Defining Wg = \og{Ig/Ic) and W\ = log(/i//c), where Ic is a normalization constant, we 
write Eq. (34) as 

riw r R 1 

(35a) 



dWg 

— cos (6*) 



ds 



l + :^7eexp(W^i) 



dWi Ai 5i 

— i = --^ + -^4exp(W^„). (356) 



We numerically solve Eq. (35) using the plasma density Ng = 3 x 10^ cm ^, the electron 
temperature T = 4200 K, frequency of the incident radiation Vg = 22.235 GHz and Ic = Ag/Bg. 
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The incident and the scattered wave intensities as functions of s, at different values of the 
scattering angle 6 = 120°, 140°, 160° and 180° in a thick plasma medium, are plotted in Fig. 2. 
At s = 0, we started with initial conditions Iq/Ic = 1 and Ii/Ic = lO"^'^. As s increases 
scattered wave intensity grows much faster than incident wave intensity, indicating transfer of 
energy from the incident wave to the scattered wave. At large s, say 3.6, and 6 = 180°, intensity 
of the incident wave becomes vary large consequently the scattered mode also becomes much 
stronger. The range of s, over which the incident and the scattered waves interact, increases 
with decreasing 9. Similarly, Fig. 3 shows the incident wave and the scattered wave intensities 
as functions of s, at different values of the scattering angle 9 = 120°, 140°, 160° and 180° in the 
case of thin plasma medium. Again at s = 0, we started with initial conditions Iq/Ic = 10^ and 
Ii/Ic = lO"^'^. At large s, say 1, and 9 = 180° intensity of the incident wave becomes vary large 
consequently scattered mode also becomes much stronger. In Figs. 2 and 3, the incident wave 
(/q) comes from the right and escaping from the boundary at s = 0, and scattered wave (/i) 
comes from the left at s = 0. The boundary of plasma at the right side exists at an arbitrary 
point, for example, at s = 3-6 in Fig. 2, or s 1 in Fig. 3. 

When the stimulated absorption of Ig and emission of Ii are very large compared to free-free 
absorption rates we can ignore the collisional damping terms in Eq. (33), therefore, we have 



Solutions of Eq. (36) can easily be obtained using the following convenient transformations-*^^ : 



- = -cos{9)-hI. 



Ol 



(36a) 




(366) 




Ux = cos{9)-^Ii. 



o 



(37) 



Now, Eq. (36) becomes 



dUo 



(38a) 



ds 
dUi 



O-l 



ds 



= UoUi. 



(386) 



In the case of forward Raman scattering, we can define a constant of motion 



Uo{s) + Ui{s) = m = Uo{0) + Ui{0) 



(39) 



the solutions of Eq. (38) can then be written as 



^■'^'- t/i(0) + S)"eip(-,».) - C""" 
Using Uo{0) = 1 and Ui{0) = 10" 4, the s dependence of the solutions is depicted in Fig. 4. We 
notice that the energy from the pump will get monotonically transferred to the scattered mode. 

B. Production of the anti-Stokes component by intersecting masers 

It is clear from the previous section that the downshifted component can be easily created 
in astrophysical conditions. However, the upshifted component can only be produced in the 
forward Raman scattering, which requires brightness temperature of the order of 10^^ K, and it 
is impossible to reach in astrophysical masers. 

We propose a special mechanism^ by invoking an inhomogeneity in the plasma to explain the 
upshifted features seen in the spectrum of NGC 4258. The possibility of enhancing the forward 
scattering by the backward scattering in the plasma, which has parabolic density profile ^^'^^ or 
by coupling Raman scattering with Brillouin scattering in the backward scattering has been 
known from the experiments on laser-plasma interactions in the laboratory plasma. 

Suppose two maser beams are intersecting in the plasma. If the maser beams are intersecting 
in such a way that they produce the same longitudinal wave at the same direction, then the 
energy transfer can take place between the two maser beams. In this process, the longitudinal 
wave produced by the backward scattering of the incident radiation is exactly frequency matched 
with the upshifted scattered wave produced by the other incident radiation. The threshold 
condition for the backward Raman scattering is 

^^thr = — l^ol'cOS^V'l > 1- (41) 
C1C3 

The frequencies of the two intersecting maser beams need not be exactly same. The threshold 
condition implies that the intersecting angle of two masers has an allowance of about ±10° 

because ki and have a beam angle of this order of magnitude. 

In this case, the transfer equations become very similar to Eqs. (24)-(26) as 

~W ^ ""d^ " ^ ^^^^ 

^ + = -C2E', - aE'^El cos (43) 

and 

^ + ^g-^ = -C3^3 + a^{EoEl cos + E'^E'l cos V^), (44) 
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where the primed quantities are for the second maser beam, and fcg — = ^3- ^^^^ process, the 
perfect phase matching is achieved for both the upshifted and downshifted components, and we 
can repeat the analysis similar to the former treatment (Sect. III). The growth rate is obtained 
as 

7 = cos ipi — COS V2 Co. (45) 

C3 \ ci ci J 

For = = 0, the threshold condition is 

{\Eo?-\E'f)^c^>l. (46) 
C3 

Therefore in order to excite the upshifted component the intensity of the first incident maser 
must be greater than the intensity of the second maser. 

IV. RESTRICTIONS OF THE THEORY 

We have not included the incoherency effect due to the incident maser radiation. The bandwidth 
and maser beam solid angle determine the coherence scale length. If we take maser beam 
solid angle to be 10~^ and the bandwidth of 1 km/sec, we can evaluate the coherency length 
approximately as Lc = A(AO)~-^ or Lc = Xc/Avo ~ 10^-10^ A, where A is the wavelength of 
the maser radiation. The Landau damping time of the longitudinal plasma wave is one order of 
magnitude higher than the coherence time scale (~ 10~^ s) of the incident radiation. Hence the 
monochromatic approximation used in the present analysis seems to be a good approximation. 
The effect of incoherence on the growth rate of stimulated Raman scattering is discussed by 
The second order Raman scattering of side band modes and the intensity level of the incident 
radiation are the main factors in limiting (saturating) the amplitude of the side band modes. 

V. CONCLUSION 

The threshold condition for backward stimulated Raman scattering is satisfied in the case of 
NGC 4258 ^. If the compact HII regions are located near the center of maser cloud, maser beam 
intersections can occur in it and the upshifted component can also be produced when the maser 
beams interact by coupling with a longitudinal plasma. The Figs. 2 and 3 indicate the length 
scales for the amplification of wave amplitudes in thick and thin plasma media, respectively. 
The inhomogeneity in the plasma and bandwidth in the incident radiation can increase the 
amplification length by one order of magnitude. 

The growth rate of stimulated Raman scattering is very large in the case of backward scat- 
tering than that in forward scattering. This lead to the complexity in explaining the spectrum 
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and time variability of extreme features. 

The red and blue shifted features could be independent and produced by two separate masers. 
A crucial test for the stimulated Raman maser model is finding the triplet with an equal sepa- 
ration and not finding a coherent time variation for the upshifted and downshifted components. 
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Figure Captions 

Fig. 1 The scattering geometry of the incident maser radiation in the case of thick plasma (a) 

and thin plasma (b) media. 
Fig. 2 The incident wave intensity lo/Ic and scattered wave intensity Ii/Ic vs s in the case 

thick plasma medium at different scattering angles = 120°, 140°, 160° and 180°. Ic is a 

normalization constant and the parameter s = {2ucUJpe/ cWo)xi. 
Fig. 3 The incident wave intensity Iq/Ic and scattered wave intensity h/Ic vs s in the case 

thin plasma medium at different scattering angles = 120°, 140°, 160° and 180°. Ic is a 

normalization constant and the parameter s = {^VpU^d clOo)x\. 
Fig. 4 Energy transfer between the coupled modes {Uq corresponding to the pump wave) and 

the parameter s = {2i'ci^pe/cuJo)xi. 
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